FIRST ORDER GLOBAL ASYMPTOTICS FOR 
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DJALIL CHAFAi, NATHAEL GOZLAN, AND PIERRE-ANDRE ZITT 

Abstract. We study a physical system of A'^ interacting particles in ffi"* subject to pair 
repulsion and confined by an external field. We establish a large deviations principle 
for their empirical measure as A'' tends to infinity. In the case of Riesz interaction, 
including Coulomb interaction in arbitrary dimension d > 2, the rate function is strictly 
convex and admits a unique minimum, the equilibrium measure, characterized via its 
potential. It follows that almost surely, the empirical distribution of the particles tends 

" ] to this equilibrium measure as A'^ tends to infinity. In the more specific case of Coulomb 

— interaction in dimension d > 2, and when the external field is a convex or increasing 

irj function of the radius, then the equilibrium measure is supported in a ring. With a 

quadratic external field, the equilibrium measure is uniform on a ball. 
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Csl 1. Introduction 

I""! We study in this work a system of N particles xi, . . . , X]\j lying in W^, d >2, for which 

Pm the energy Hj^{xi, . . . , xn) of the configuration (xi, . . . , xn) G (M'^)^ takes the form 
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Hn{xu ...,xn):= NY^Vixi) + ^ W{xi, Xj) 

i=l i<j 



where 1/ : M'^ — ^ M and 1/^ : M^ x M"^ — ^ (— oo, -|-cxd] is finite outside the diagonal, and sym- 
metric: W{x,y) = W{y,x) for all x,y £ M.'^. We can express the energy Hn{xi, . . . ,xn) 
j> of the configuration in terms of the empirical distribution of the particles 

v^ 1 ^ 



N 



■^ Namely, we simply have 

^ Hn{xi,...,xn) = N'^ (p{x,y)diJ.N{x)dfiN{y) + N V{x)dfiN{x) 

• • '7~y 

> 

• ^H where the double integral is off-diagonal and where 99 : M x M — ?► R U {+00} is given by 

^ ^{x,y):=^{V{x) + V{y)+W{x,y)). 

Our work is motivated by the following physical control problem: given the "interaction 
potential" W, for instance a Coulomb potential, and a target probability measure ^^ on 
M^, for instance the uniform law on the unit ball, can we tune the "external localizing 
potential" (or "external localizing field") V such that hn —>■//* as A'' —)■ 00? In this 
direction, we provide some partial answers in Theorem 1.1, Theorem 1.2 and Example 1.3 
below, and also discuss some possible extensions. 
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1.1. Main results. From now on, and unless otherwise stated, we denote by |-| the Eu- 
clidean norm of M'^ and we make the following additional assumptions: 

(HI) The function ty : M'^ x R'^ — ;► Ruj+oo} is continuous on R^ x M.'^, symmetric, takes 
finite values on R x M \ {(x, x); x G R } and satisfies the following integrability 
condition: for all compact set K cM. , the function 

axiu) = sup{W{x, y);\x-y\ >u,x,y£ K}, 

with the convention ax^u) = if u > diam.(K), is such that the function x M- 
ai^(|x|) is locally Lebesgue-integrable on R*^. 
(H2) The function y : R — ^ R is continuous and such that lim|2,|_i.^oo V{x) = +00 and 



/ exp {—V{x)) dx < 00. 



(H3) For some £« e (0, 1/2), the function {x,y) ^ (l - ^o) iV{x) + V{y)) + lW{x,y) 
is bounded from below. 

It is easily seen that under (H2)-(H3), for ah N £ W , 

Zn ■= ■ ■ ■ / 6xp {—H]\[{xi, . . . , X]\f)) dxi ■ ■ ■ dx]\f < cxd, 

JR'i jRd 

so that we can define the Boltzmann-Gibbs probability measure P/v on (R'^)^ by 

dPN{xi,...,XN) := -^e-"^^'''^-'^^Uxi...dxN. (1.1) 

This is the distribution of a system of A^ particles living on R^, subject to a confining 
potential NV acting on each single particle and to an interaction potential W acting on 
each pair of particles. Following [ForlO, Section 11.3.1] (see also [DXOl, Section 9.6]), the 
model is of Calogero-Sutherland type, and the density of P/v takes the form 

N 

Y[fi{xi) n h{xi,Xj). 

i=l l<i<j<N 

Let A^i(R ) be the set of probability measures on R . The mean-field symmetries of the 
model suggest to study, under the exchangeable measure P^, the behavior as A^ — ;► cxd of 
the empirical measure ^^v, which is a random variable on A^i(R ). With this asymptotic 
analysis in mind, we introduce the functional / : A^i(R ) — t- (—00, -|-oo] given by 

lifJ-) ■= hp{x,y)dfj.{x)dn{y). 

Note that (H2) and (H3) imply that ip is bounded from below, so that this integral makes 
sense (see Lemma 2.2). 

Theorem 1.1 below is our first main result. It is of topological nature, inspired from 
the available results for logarithmic Coulomb gases in random matrix theory [P)AG97, 
BAZ98, PII98, Harl2]. We equip Mi(R'^) with the weak topology (defined by duaHty 
with bounded continuous functions), and for any set A C A^i(R'^) we denote by int(^), 
clo(^) the interior and closure of A with respect to this topology. 

To formulate the large deviation result we need to introduce the following additional 
technical assumption: 

(H4) For all u € 7Wi(R'^) such that I{u) < -|-oo, there is a sequence {v'n)nm of prob- 
ability measures, absolutely continuous with respect to Lebesgue, such that i/„ 
converges weakly to z/ and I^i'n) — ^ -^(^)) when n — )• 00. 
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It turns out that assumption (H4) is satisfied for a large class of potentials V, W, and 
several examples are given in Proposition 2.8 and Theorem 1.2. We recall that the weak 
topology on A4i(JR.'^) can be metrized [Dud02] by the Fortet-Mourier distance defined by 



difi,i^):= sup \ fdfi- fdi^\. (1.2) 

max(|/|„o,|/|L.p)<l U J J 

Theorem 1.1 (Large Deviations Principle). // (H1)-(H2)-(H3) are satisfied then 

(1) I is lower semi- continuous and has compact level sets; 

(2) Under {Pn)n, the sequence {^j.n)n of random elements of Aii{W^) equipped with 
the weak topology has the following asymptotic properties. For every Borel subset 
A ofMiiW'), 

limsup— ^logZArPAr(//Ar ^ A) < - inf I{^) 

and 
— inU I (fi); fj, E mt( A), fi <^ Lebesgue} < liminf — ^\og Z^P^I^n G A). 

(3) Under the additional assumption (H4), the full Large Deviation Principle at speed 
I/N"^ holds with the rate function 

h:=I - inf /. 

Xi(M'*) 

More precisely, for all Borel set A C 7Wi(M ), 

- inf h{fi) < lim inf —^ log Pn{hn & A) 

< lim sup -— T log P^ifJ-N ^ A) < - inf I^(/i). 

Af^oo A* /ieclo(A) 

(4) Let Imin '■= {^ S A4i : /*(/u) = 0} 7^ 0. // (H4) is satisfied and if {iJ'n)]\[>i o-re 
constructed on the same probability space, and if d stands for the Fortet-Mourier 
distance (1.2), then we have, almost surely, 

lim d{fJ,N,Imin) = 0. 

Af— s>co 

A careful reading of the proof of Theorem 1.2 indicates that if Imin = {/^*} is a singleton, 
and if (H4) holds for i/ = //^, then ^^v -^ /^* almost surely as A'' — >■ oo. 

Our second main result, expressed in Theorem 1.2 and Example 1.3 below is of dif- 
ferential nature. It is based on an instance of the general Gauss problem in potential 
theory [Zor03a, Zor03b]. It concerns special choices of V and W for which I^ achieves its 
minimum for a unique and explicit //^ S A4i(M ). Recall that the Coulomb interactions 
correspond to the choice W{x,y) = k/\(x — y) where A;^ is the Coulomb kernel (opposite 
in sign to the Newton kernel) defined on M , d> 1, by 

■\x\ if d = 1, 
feA(x):=<|log^ if(i = 2, (1.3) 

This is, up to a multiplicative constant, the fundamental solution of the Laplace equation. 
In other words, denoting A := 9^^ + • • • + d"^ the Laplacian, we have, in a weak sense, in 



There is no boundary conditions here, and thus the term "Green function" is not appropriate. 



4 DJALIL CHAFAI, NATHAEL GOZLAN, AND PIERRE-ANDRE ZITT 

the space of Schwartz-Sobolev distributions D'(M'^), 

fi ifd = l, 

- cA/cA = So with c := < ^ if d = 2, (1.4) 

,, , „s — if d > 3, 

where w^ := p^ ^,^-^ is the volume of the unit ball of W^. Our notation is motivated by the 
fact that —A is a nonnegative operator. The case of Coulomb interactions in dimension 
d = 2 is known as "logarithmic potential with external field" and is widely studied in the 
literature: see [HPOO, ST97, AGZIO, Harl2]. To focus on novelty, we will not study the 
Coulomb kernel for d < 2. We refer to [Len63, EL62, Len61, BL75, AM81, SS13] and 
references therein for the Coulomb case in dimension d = 1, to [BAG97, AGZIO, Harl2] 
to the Coulomb case in dimension d = 2 with support restriction on a line, to [BAZ98, 
PH98, HPOO, Harl2, ST97, SS12] for the Coulomb case in dimension d = 2. 

The asymptotic analysis of /Utv as A^ — )• oo for Coulomb interactions in dimension 
d > 3 motivates our next result, which is stated for the more general Riesz interactions in 
dimension d > 1. The Riesz interactions correspond to the choice W{x,y) = k^^{x — y) 
where /^Aqi < a < d, d > 1, is the Riesz kernel defined on W^, by 

k^^{x):=-^. (1.5) 

\x\ 

Up to a multiplicative constant, this is the fundamental solution of a fractional Laplace 
equation (which is the true Laplace equation (1.4) when a = 2), namely 

- cA„/cA„ = ^-^1) = <5o with c:=^^^^, (1.6) 

where the Fourier transform T and the fractional Laplacian Aq, are given by 

^(A;aJ(0:= / e'*"^"A:Ajx)dx = — 4^ and ^o.f := -^'K^J'~\\£.rHf))- 

Note that A2 = A while Aq, is a non-local integro-differential operator when a 7^ 2. 
When d > 3 and a = 2 then Riesz interactions coincide with Coulomb interactions and 
the constants match. Beware that our notations differ slightly from the ones of [Lan72, 
p. 44] . Several aspects of the Gauss problem in the Riesz case are studied in [DS07, Zor03a, 
Zor03b]. 

In the Riesz case, to any reasonable probability measure /i, corresponds a function U^ 
called the potential of /i — its precise definition is recalled below, see Section 3. One can 
also define a notion of capacity of sets, and a property is said to hold quasi everywhere if 
it holds outside a set of zero capacity (once more, see Section 3 for details). The following 
theorem is essentially the analogue in M of a result of Dragnev and Saff on the sphere 
(see [DS07]). 

Theorem 1.2 (Riesz gases). Suppose that W is the Riesz kernel W{x, y) = (3k^^{x — y), 
for some (3 > 0. Then: 

(1) The functional I is strictly convex; 

(2) (HI )- (H2)- (H3)- (H4) are satisfied and Theorem 1.1 applies; 

(3) There exists a unique fii, G MiiW^) such that 

Hi^*)= inf ^,^(/^); 



(4) If we define {fJ-N)N on a unique probability space, then with probability one, 

hm UN = fJ-^- 
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If we denote by C^, the real number 

C, = J (C/^^ + V)dfi, = J(/i,) + J Vdfi,, 

then the following additional properties hold: 

(5) The minimizer ^^ has compact support, and satisfies 

U^*{x) + V{x) > Ci, quasi everywhere, (1-7) 

U^'*{x) + V{x) = C^ for all X £ supp fi. (1.8) 

(6) If a compactly supported measure /i creates a potential U^ such that, for some 
constant C S M, 

U^{x) + V{x) = C on supp(/i), (1.9) 

f/M _|_ y > (^ quasi everywhere, (1-10) 

then C = Ci, and fi = fj.-),. 

The constant C* is called the "modified Robin constant", see e.g. [ST97], where the 
properties (1.7-1.8) and the characterization (1.9-1.10) are established for the logarithmic 
potential in dimension 2. The minimizer fj,^ is called the equilibrium measure. 

Example 1.3 (Equilibrium measure of Coulomb gases with radial external fields). Suppose 

that for a fixed real parameter /3 > 0, and for every x, y G M , 

V{x) = v{\x\) and W{x,y) = /3k/\(x — y), 

where v is two times differentiable. Denote by da^- the Lebesgue measure on the sphere of 
radius r, and let ad be the total mass of dai. Let w{r) = r'^~^v'{r), and suppose either 
that V is convex, or that w is increasing. Define two radii rg < Rq by: 

ro = inf {r > 0; v'{r) > 0} and w{Rq) = p{d - 2) . 

Then the equilibrium measure fii, is supported on the ring {x; \x\ G [ro,i?o]}, o^nd is abso- 
lutely continuous with respect to Lebesgue measure: 

w' (r) 
dfj,{r) = M{r)dardr where M{r) = .^ _ ^, ^^_^ \^^^R^^^{r). 

In particular, when v{t) = t^ then ^^ is the uniform distribution on the centered ball of 
radius (^(d - 2)/2)V'i. 

This result on Coulomb gases with radial external fields belongs to the folklore of Poten- 
tial Theory, and can be found for instance in [LGIO, Proposition 2.13]. It follows quickly 
from the Gauss averaging principle and the characterization (1.9-1.10). For the sake of 
completeness, we give a (short) proof in Section 4-3. 

In particular, by using Theorem 1.2 with a = 2 together with Example 1.3, we obtain 
that the empirical measure of a Coulomb gas with quadratic external field tends almost 
surely to the uniform distribution on a ball when N ^ oo. This phenomenon is the 
analogue in arbitrary dimensions of the well known result in dimension 2 for the logarithmic 
potential with external field (Complex Ginibre Ensemble, [BAZ98, PH98]). 

1.2. Outline of the article. In the remainder of this introduction, we give several com- 
ments on our results, their links with different domains, and possible directions for further 
research. The large deviations principle stated in Theorem 1.1 is proved in Section 2. 
The proofs of Theorem 1.2 and Example 1.3 use several concepts and tools from Potential 
Theory, which we recall synthetically and discuss in Section 3 for the sake of clarity and 
completeness. The proofs themselves are in the final section. 

1.3. Comments, possible extensions and related topics. 
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Equilibrium measure for Riesz interaction with radial external field. To the knowledge of 
the authors, the computation of the equihbrium measure for Riesz interactions with radial 
external field, beyond the more specific Coulomb case of Example 1.3, i.e. for a 7^ 2, is an 
open problem, due to the lack of the Gauss averaging principle when a ^ 2. 

Beyond the Riesz and Coulomb interactions. Theorem 1.2 concerns the minimization of the 
Riesz interaction potential with an external field V, and includes the Coulomb interaction 
if d > 3. In classical Physics, the problem of minimization of the Coulomb interaction 
energy with an external field is known as the Gauss variational problem [Zor03a, Zor03b]. 
Beyond the Riesz and Coulomb potentials, the driving structural idea behind Theorem 1.2 
is that if W is of the form W{x,y) = kL{x — y) where k^ is the fundamental solution 
of an equation dt — L = where L is a local (say differential) operator, and if V is 
super-harmonic for L, i.e. LV > 0, then the density of ^u* is roughly given by LV up to 
support constraints. This can be easily understood formally with Lagrange multipliers. 
The limiting measure //^ depends on V and W, and is thus non-universal in general. In 
the opposite direction, if W and ^u* are prescribed, and if // M- ffW{x,y) d^{x)dfj,{y) is 
convex, then the proof of Theorem 1.2 suggests a natural candidate for the external field: 
V{x) = —Uy^{x) := Ci, — fW{x,y) dfi^,{y) for all x G supp(/i^), where C* is a constant. 

Second order asymptotic analysis. The asymptotic analysis of fiN — /U^ as A^ — )■ 00 is a 
natural problem, which can be studied on various classes of tests functions. It is well 
known that a repulsive interaction may affect dramatically the speed of convergence, and 
make it dependent over the regularity of the test function. In another direction, one may 
incorporate an inverse temperature parameter /3 in front of the energy H]\f, and study 
the low temperature regime, i.e. asymptotics as /3 — >■ 00. In the Coulomb case, this leads 
to Fekete points. The analysis of the second order when both /3 — )■ 00 and A^ — ;► cxd is 
conjectured to be related to an Abrikozov lattice, see [SS12, SS13], and [Bctl3]. 

Edge behavior. Suppose that V is radially symmetric and that fi^, is supported in the 
centered ball of radius r, like in Example 1.3. Then one may ask if the radius of the 
particle system maxi<fc<„ jx^l converges to the edge r of the limiting support as A^ — >■ 00. 
This is not provided by the weak convergence of fi^- The next question is the fluctuation. 

Topology. It is known that the weak topology can be upgraded to a Wasserstein topology 
in the classical Sanov theorem for empirical measures of i.i.d. sequences, see [WWWIO], 
provided that tails are strong exponentially integrable. It is then quite natural to ask 
about such an upgrade for Theorem 1.1. 

Connection to random matrices. Our initial inspiration came, when writing the survey 
[BC12], from the role played by the logarithmic potential in the analysis of the Ginibre 
ensemble. When d = 2 and V(x) = \x\ and W(x,y) = I3k^(x — y) = /31og , j_ , with 
(3 = 2 then Pjv is the law of the (complex) eigenvalues of the complex Ginibre ensemble: 

dPN{x) = Z^^e-^S-i l^'l' Jl \xi - Xj\^dx. 

i<j 

(here M^ = C and P/v is the law of the eigenvalues of a random N x N matrix with i.i.d. 
complex Gaussian entries of covariance ^It)- For a non-quadratic V , we may see P^ as 
the law of the spectrum of random normal matrices such as the ones studied in [AHMll]. 
On the other hand, in the case where d = 1 and Vix) = Ixl and Wix,y) = /31og 1 _ 1 
with /3 > then 

dPN{x) = ZJ^e-^^ti l^^l' [] \xi - xjf dx. 

i<j 

This is known as the /3-Ensemble in Random Matrix Theory. For /3 = 1, we recover the 
law of the eigenvalues of the Gaussian Orthogonal Ensemble (GOE) of random symmetric 
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matrices, while for /3 = 2, we recover the law of the eigenvalues of the Gaussian Unitary 
Ensemble (GUE) of random Hermitian matrices. It is worthwhile to notice that —log |-| 
is the Coulomb potential in dimension d = 2, and not in dimension d = 1. For this reason, 
we may interpret the eigenvalues of GOE/GUE as being a system of charged particles in 
dimension d = 2, experiencing Coulomb repulsion and an external quadratic field, but 
constrained to stay on the real axis. We believe this type of support constraint can be 
incorporated in our initial model, at the price of a bit heavier notations and analysis. 

Non- compact equilibrium m,easures. The assumptions made on the external field V in 
theorems 1.1 and 1.2 explain why the equilibrium measure ;U^ is compactly supported. If 
one allows a weaker behavior of V at infinity, then one may produce equilibrium measures 
//j, which are not compactly supported (and even heavy tailed probability measures). This 
requires to adapt some of the arguments, and one may use compactification as in [liar 12]. 

Sim.ulation problem.. It is natural to ask about the best way to simulate the probability 
measure Pn. A pure rejection algorithm is too naive. Some exact algorithms are available 
in the determinantal case d = 2 and W{x, y) = —2 log |x — y| , see [HKPV06, Algorithm 18] 
and [SZT09] . One may prefer to use a non exact algorithm such as a Hastings-Metropolis 
algorithm. One may also use an Euler scheme to simulate a stochastic process for which 
Pm is invariant, or use a Metropolis adjusted Langevin approach (MALA) [RROl]. 

More general energies. One may study more general energies with many bodies interac- 
tions, of the form, for some prescribed symmetric W^ : (JR.'^)^ i— t-M, 1 < k < K, K > 1, 

K 

Fjv(xi,...,xjv) = ^iV2-fc J2 Wkixi,,...,XiJ. 

k=l ii<--<ik 

This leads to the following candidate for the asymptotic first order global energy functional: 

/iH> ^2"^ /••• W{xi, . . . ,Xk) dfi{xi) ■ ■ ■ dfi{xk). 
k=i -^ -^ 

Stochastic processes. Under general assumptions on V and W ., see for instance [Roy07], 
the law Pat is the steady state solution of the Fokker-Planck evolution equation dt — L = 
where L is the linear elliptic Markov diffusion operator 

L := A - VHn • V, 

acting as Lf = A/— (Vf^A? , V/) . The operator L is self-adjoint in L^(P/v). In probabilistic 
terms, the law Pjv is a reversible invariant probability measure for the Markov diffusion 
process (X()j>q solution of the system of Kolmogorov stochastic differential equations 

dXt = V2dBt- VHN{Xt)dt 

where (i?i)^>Q is a standard Brownian motion on (M ) . In the case where V{x) = \x\ 
and W = (no interaction) then P/v is a standard Gaussian law M{0,IdN) on (M'^)^ 
and (Xt)j>Q is an Ornstein-Uhlenbeck Gaussian process; while in the case where d = 1 

and V{x) = \x\ and W{x,y) = —f3log\x — y\ then P/v is the law of the spectrum of 
a /3-Ensemble of random matrices and (X()^>q is a so called Dyson Brownian motion 
[AGZIO]. If /XAT^i is the law of Xt then ^.n^yf^t — ^ IE/xtv weakly as t — )■ cxd. The study 
of the dynamic aspects is an interesting problem connected to McKean-Vlasov models 
[CL97, Fon04, LLX13]. 
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Calogero-(Moser-) Sutherland Schrodinger operators. Let us keep the notations used above, 
and let us consider the isometry Q : L^(P/v) — )■ L'^{dx) defined by 



e(/)(x) := f{x)J^^^ = /(:c)e-i(^-(-)+i°s(^-)). 

The differential operator S := —QLQ^^ is a Schrodinger operator: 

S:=-eLQ-^ = -A + U, U :=hvHNf -^AHn 

which acts as Sf = —A/ + Uf. The operator S is self-adjoint in L'^{dx). Being isometri- 
cally conjugated, the operators —L and S have the same spectrum, and their eigenspaces 
are isometric. In the case where V{x) = \x\ and W = (no interactions), we find that 
and U = 2(1 — V), and S is a harmonic oscillator. On the other hand, following [ForlO, 
Proposition 11.3.1], in the case d = 1 and W{x,y) = —log \x — y\ (Coulomb interaction), 
then S" is a Calogero-(Moser-)Sutherland Schrodinger operator: 

^- A F+^Vr2 ^ V 1 _iV N{N-1) 
S--A-Eo + -}_^x,-^ Z. (x-x-)^' ^°-y + 2 ■ 

^ i=l ^ l<i<j<N ^^* ^^> ^ ^ 

More examples are given in [ForlO, Proposition 11.3.2], related to classical ensembles of 
random matrices. The study of the spectrum and eigenfunctions of such operators is a 
wide subject, connected to Dunkl operators. These models attracted some attention due 
to the fact that for several natural choices of the potentials V^ W, they are exactly solvable 
(or integrable). We refer to [ForlO, Section 11.3.1], [DXOl, Section 9.6], [CDR+08, Section 
2.7] and references therein. 

2. Proof of the large deviations principle — Theorem 1.1 
The proof of Theorem 1.1 is split is several steps. 

2.1. A standard reduction. To prove Theorem 1.1, we will use the following standard 
reduction explained in the following proposition (see for instance [DZ9N, Chapter 4]). 

Proposition 2.1. Let {Qn)n£N* be a sequence of probability measures on some Polish 
space {X,d), (Ziv)^gpj. and {eN)NGN* two sequences of positive numbers with en ^ and 
X : A' — 7- M be a function bounded from below. 

(1) Suppose that the sequence {Qn)n£W satisfies the following conditions 

(a) The sequence (ZjvQjv)AfeN* is exponentially tight: for all L >0 there exists a 
compact set Kl C X such that 

lim sup En log Z]\fQN {X \ Ki) < —L. 

Af-s-oo 

(b) For all x e X, 

lim limsupffAT log Ziv(5Ar(-B(x, r)) < —I{x). 

Then the sequence {ZiyQj^)j\f^fq* satisfies the following large deviation upper bound: 
for all Borel set A C X , it holds 

limsupeNlog ZnQn{A) < -inf {X(/i;/i e clo(^))}. (2.1) 

(2) If in addition [ZnQn)^^^* satisfies the following large deviation lower bound: for 
all Borel set Ad X, 

-inf{X(x);a; G int(^)} < liminf e7vlogZArQ7v(A), (2.2) 
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then {QN)NeN* satisfies the full Large Deviation Principle with speed em and rate 
function X^ = X — m.ix(^x^{x), namely for all Borel set A d X , 

— inf{X^(x);x G int(j4)} < liminf eArlog(5Ar(j4) 

A''— s>oo 
< limsupeArlogQAr(yl) < — inf {X-^(3;); x G c\o[A)}. 

Af-5>oo 

Proof. Let us begin by (1). Let (5 > 0; by assumption, for any x G Af, there is jy^; > such 
that 

hmsupeArlogZAr(57v(5(a;,??a;)) < -"^{x) +5. 

If F C A" is compact, there is a finite family (xi)i<i<m of points of F such that F C 
U^iB(xi,7?^J. Therefore, 

hmsupeArlogZArQAr(F) < hmsup^Arlog ^ZNQN{B{xi,rixJ) 

= max limsupENlog {ZNQN{B{xi,i]xi))) 

l<i<m 7V->oo 

< max — X(xj) + 6 < — inf X + 6. 

l<i<m F 

Letting 5—7-0 yields to (2.1) for A = F compact. 

Now if F is an arbitrary closed set, then for all L > 0, since FdICl is compact, it holds 

lim sup En log Z]\[Qn{F) 

< max f limsupeArlogZAfQivC^ H Kl), limsup EN^og Z]\[Q]\j{Kl) 

< maxf— inf X\—L). 

FHKl 

Letting L — t- oo shows that (2.1) is true for arbitrary closed sets F. Since A C clo(A), the 
upper bound (2.1) holds for arbitrary Borel sets A. 
To prove (2), take A = X in (2.2) and (2.1) to get: 

lim En log(ZAr) = - inf X G M. 
Subtracting this to (2.2) and (2.1) gives the LDP with rate function X^. D 

In our context, X = 7Wi(M ) is equipped with the Fortet-Mourier distance (1.2). 

2.2. Properties of the rate function. Let us prove the properties of the rate function 
/j. announced in Theorem 1.1, point (1). 

Lemma 2.2 (/^ is a good rate function). Under Assumptions (H2)-(H3), the functional 
/j, is well defined and is a good rate function (the levels sets {I < k} are compact). 

Proof. According to Assumption (H3), there exists Eo G (0, 1/2) such that 

^{x,y)>Eo{V{x) + V{y)) + c, (2.3) 

with c = inf2:^y(l/2 — Eo)iy{x) + V{y)) + {l/2)W{x, y) > — oo. It follows from Assumption 
(H2) that if is bounded from below. The functional I* is thus well defined with values 
in [0,oo]. Furthermore, according to the monotone convergence theorem, / = sup„gpj/„, 
where In{f^) = JJ min(i^(x, y); n) d^{x)djjL{y). The functions min(c/9, n) being bounded and 
continuous, it follows that the functionals /„ are continuous for the weak topology (see for 
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instance [DZ98, Lemma 7.3.12]). Being a supremum of continuous functions, / is lower 
semi-continuous. Set 6^ = inf (/?; we have, for every fj, G A^i(M"'), L > 0, 



I{n) -K= {<f{x, y) - K) dn{x) dfi{y) 

> / / '^\x\>L,\y\>L{v{x, y) - K) dfi{x) dfi{y) 

>{bL-K)K\x\ >Lf. 

where hi '■= inf\ri.\>L.\y\>L^ix,y). According to (2.3) and Assumption (H2), we see that 
bL — ^ +00 as L — ^ +00. Therefore, there exists L^ > such that 6/, > 6^ for every L > L^. 
We get then for every real number r >h^, 

{^i G MiiM.'^) : I{p) < r} C |/x G Mi{W') : ^i{\x\ > L) < J^^^,L > lA. 

Since bi — >■ +00 as L — >■ +00, the subset of M-i{W^) in the right hand side is tight, 
and the Prohorov theorem implies then that it is relatively compact for the topology of 
A4i(R'^). Since / is lower semi-continuous, the set {/ < r} is also closed, which completes 
the proof. D 

2.3. Proof of the upper bound. For all A^ G N*, one denotes by Qn the law of ^n = 
jf J2i=i ^Xi under the probability P/v defined by (1.1): Qn is an element of A^i(A^i(M'^)). 

Lemma 2.3 (Exponential tightness). Under Assumptions (H2)-(H3), the sequence of 
measures {ZNQN)Nm* is exponentially tight: for all L > there exists a compact set 
Kl C Mi{W^) such that 

limsupeTvlogZArQTv {Mi{M.'^) \ kA < -L. 

Proof Let us consider the probability uy G A^i(M ) by vv{dx) = tjC" '^' dx, with C = 
J e-^(^) dx. For all L > 0, 



ZnQn i V dfiN > L 



= C'j\fy .,.>.} exp [-N^ Jl^^ ^ d,%) d4 

<C'' J l{j v'dM^>L} exp (-A2 JJ^ eo{V{x) + Viy)) + 6(i/x^(x)(i/x^(y)) du^ 

= C''e-'>^(^-^^Jl^jy,^^^,^e.v {-2N\l - N--)e. j V d^.^) d.^ 

where the first inequality follows from (2.3), with Eq G (0, 1/2) and 6 G M independent of 
N . This implies that 

limsup — -T log ZnQn [ / V dfiN > L] <b- 2eoL. 

Since V{x) — ?► cx) when |x| — t- -|-oo and V is continuous, it follows from Prohorov theorem 
that the sets K^ := {/x G Mi{M.'^); J V dfj, < L}, L > are compact. This completes the 
proof. n 

Proposition 2.4 (Upper bound). Under Assumptions (H2)-(H3), for all r > 0, for all 

11 G Mi{W^), 

lini limsup—^ log ZAfQAr(S(/i,r)) < -/(/i). 
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where the ball B(fj,,r) is defined for the Fortet-Mourier distance (1-2). 

Proof. Recall that vv{dx) = C'^e-^^^) dx and let /x e Mi{W^), r > and n G N*; it 
holds 

ZNQN{B{p,,r)) 

= C^ li?(^,r)(AtAf)exp f-iV^ // ip(x,y) d^iN{x)diiN{y)\dvy 

<C^ j lB(^,r)(AtAf)exp f-iV^ // min((/?(x,2/);n)d^Af(a;)d^Af(y) 

<C^j lB(^,r)(AtAf)exp f-iV^ // m.m.{Lp{x,y)]n) diiN{x)diiN{y) 

X exp A^ / min((^(x,x); n) (i/iAr(x) Idz^y 



dv^ 



< C^e"^ exp -iV^ jj^f / / min(c^(x, y); n)(ii/(x)dz^(y) . 

\^ ueB{fi,r)Jj J 

Therefore, 

limsup— ^ log ZNQN{B{fj,, r)) <- inf /„(j^), 

where /n(j^) := JJ min(ip{x,y);n)di'{x)du{y). Arguing as in the proof of Lemma 2.2, we 
see that /„ is bounded continuous. Therefore, limr_5.oinf,ygB(^^r) -^n(z^) = InifJ')- Finally, 
the monotone convergence theorem implies that sup„gpj. InifJ-) = -^(/u), which ends the 
proof. □ 

Using this Proposition, Lemma 2.3 and the first point of Proposition 2.1, we get the 
upper bound of Theorem 1.1, point (2). 

2.4. The lower bound and the full LDP. In what follows, we denote by \A\ the 
Lebesgue measure of a Borel set A C M". 

Proposition 2.5 (Lower bound for regular probabilities). Under the assumptions (Hl)- 
(H2)-(H3), for all probability measure jj, on W^ supported in a box B = ni=i['^i) ^«]; 
Oj, 6j G M, with a density h with respect to the Lebesgue measure such that, for some 5 > 0, 
5 < h < 6^^ on B it holds 

N^oo iV ^ 

where B{^, r) denotes the open ball of center ji and radius r for the Fortet-Mourier dis- 
tance (1.2). 

If B is the box Hi'^fei&fc] in ^'^^ let KB) and L{B) be the minimum (resp. maximum) 
edge length: 

1{B) = min(6fc - at), L{B) = max(6fc - a^) 

k k 

We admit for a moment the following result: 

Lemma 2.6 (Existence of nice partitions). For all d and all 5 > there exists a constant 
C{d, 5) such that the following holds. For any box B, any integer n, and any measure /i 
with a density h w.r.t. Lebesgue measure, if S < h < 6^^, then there exists a partition 
{Bi, B2 ■ ■ . Bn) of B in n sub-boxes, such that: 

(1) B is split in equal parts: for all i, id{Bi) = ^^{B) ; 

(2) The edge lengths of the Bi are controlled: 



1 u„.^„„^^,,„^^C{d.S) 



oiijW'"-''^ < im < HB.) < 



n 



l/d 



L(B). 
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Proof of Proposition 2.5. For each A'^ we apply Lemma 2.6 to obtain a partition of B in 
A^ boxes B^ , . . . B^. We let djv be the maximum diameter of the boxes: by the lemma, 
since IJ,{B) = 1, 

-j^ - ^^^^^^ '^^ - i^' 

where ci and C2 only depend on B, d and 5. 

Note that, for all 1-Lipschitz function / with \\f\\^ < 1, if Xi G -Bj- for all i < N, since 
fJ-^Bi) = 1/N we have: 

1 ^ Af 



i=l 



i=l 



Therefore, if A'^ is large enough, 

B^ x■■■xB^Cl{^iN^B{^i,r)}. 

Let us denote by C/^ C B^^ the box obtained from B^^ by an homothetic transformation 
of center the center of B^ and ratio (say) 1/2. It holds. 



N \ N 

N\ 
i 



ZNQN{B{ti,r)) > exp -iY^maxF - J2 ^^x H^ [] 1^/ 



i=l ^i i<j ^i ^S / i=l 



Since 



C 



N 



> (/(i?jY)/2) > cs/N for some absolute constant C3, we conclude that 



1 1^1 

liminf -— TT log (ZjvQAf(i?(w,r))) > — limsup -— > maxl^ — limsup —: TT > max W. 

For all A^, consider the locally constants functions V/v : -B — ?■ M and Wn : B x B ^ M 
defined by Vn{x) = max^iv V, ii x & B^ and WN{x,y) = max^^jv^^iv W, if x S B^ and 

y G Bf . Since ;u(Sf ) = l/A^, it holds 

1 ^ f 1 1 /■ 

— VmaxF=/ VN{x)d^J.{x) and ^r^o X! S^^^., ^ = o / WAr(x, y) (i^(x)(i/i(y). 

The uniform continuity ol V on B immediately implies that Vn converges uniformly to 
V, and so / Vjy dfj. ^ J V dfi. For the same reason Wn converges uniformly to W on 
{B X B) n {{x,y) G M'^ X R'^;\x-y\ > u}, for ah u > 0. To apply the dominated 
convergence theorem it remains to bound Wn from above by some integrable function. 
Recall that asiu) = sup\^_y\y^W{x,y), so that W{x,y) < obHx — y\). Obviously 

max \x — y\<2d]s[+ min Ix — yl. 

By construction, since i / j, min^^y^^^N^c'^ \x - y\ > \{l{Bf) + l{Bf)) > ^N''^/'^ > 
^djy- Therefore, the right hand side is bounded by (8C2/C1 + 1) minj^Ar^civ \x — y\, so 
there is an absolute constant C4 such that 

min \x — y\ > C4, max \x — y\. 

(x,j/)eCfxCf (x,y)eBfxi?f 

Since the function as is non-increasing, it holds 

max aB(|x-y|) < min aB{cA\x-y\). 

s that Wn{x, y) < c^OBic^lx — y\), x j^ y, for si 
C5, and using Assumption (HI), it is not difficult to see that the function in the right 



We conclude from this that Wn{x, y) < c^OBic^lx — y\), x j^ y, for some absolute constant 
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hand side is integrable on B with respect to /i x /i. Applying the dominated convergence 
theorem, we conclude that 

linrinf ^ log {Z^QNiBifi, r))) y-J V{x) dis{x) -^JJ W{x, y) df,{x)di2{y) 

Let us now prove that "nice" partitions exist. 

Proof of Lemma 2. 6. The proof is an induction on the dimension d. 

Base case. Let d = 1, and suppose that B = [ao,6o]- Since (i has a density, there 
exists "quantiles" oq = go < (Zi • • • < 9n = ^o such that 

VI < i < n, n{[qi^i,qi]) = -ii{B). 

n 

In this simple case 1{B) = L{B) = bo — oq, and l{Bi) = L{Bi) = (qi — qi-^i). The 
boundedness assumption on h implies that 

1 

5{qi - qi-i) < fi{[qi^i,qi]) < -{qi - qi-i) 

6{bo - ao) < fi{B) < t(&o - ao), 

and the claim holds for d = 1 with C(l, 5) = p . 

Induction step. Suppose that the statement holds for a dimension d — 1. Let B = 
[oq, bo] X B' he a box in dimension d (where B' is a (d — l)-dimensional box). Let /io be 
the first marginal of ^ (this is a measure on [ao, bo] C M). 

Let b = [n^''^\ be the integer part of n^''^. If b = 1, n < 2'^, and we can reason as in 
the base case, on the one-dimensional measure fio, to find a partition of i? in n slices of 
mass fj,{B)/n. Since the number of slices is less than the constant 2^^, the edge length is 
controlled as needed. 

If 6 > 2, we look for a decomposition of n as a sum of b integers rii, each as close to 
^{d-i}/d g^g possible: the idea is to cut B along the first dimension in b slices, and to apply 
the induction hypothesis to cut the slice i in rii parts. 

To this end, decompose the integer n in base b: 

d-l 

3ao, ai, . . . ad-i G {0, ... 6 - 1}'', n = b'^ + Y, (^kb^ ■ 

fc=0 

Writing at = Ei=i l{i<afc} we get 

b / d-l \ b 

1=1 \ k=0 / i=l 

where ni = b'^~^ + Yl,k=o '^{i<ak}^^ ■ From this expression, we get the bound b'^~^ ^ rn < 
{b'^ — l)/(6 — 1). Since k — 1 > k/2 whenever /c > 2, using the inequalities b < v}''^ and 
b > n^/^ — 1 > ^n^^'^, we get 

J_^{d-i)/d < < 2n('^-i)/^. 
2d-i - - 

Now let us cut B along its first dimension. Recall that ^o is the first marginal of fi. By 
continuity there exist quantiles ao = qo < qi < ■ ■ ■ qt = bo such that 

VI < i < 6, fii{[qi-i,qi]) = n{[qi-i,qi] x B') = — /x(-B). 

n 
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We apply the induction hypothesis separately for each 1 < i < 6, to the {d— l)-diniensional 
box B' , with the measure 

Mi(-) = Ai(te-i,'7i] X •) 
and the integer m to obtain a decomposition B' = lj?=i B[ • such that: 

(1) the edge lengths B[ , are controlled, 

(2) ^i,{B[^^) = lpi,{B'). 

Finally, for all 1 < i < 6 and all 1 < j < Ui, let 

Bi,j = [qi-i,qi] xB'ij. 
Let us check that the partition B = [J^ U • Bij satisfies the requirements. By definition, 

KBij) = fJ'iiB'ij) = —Hi{B') = -fi{B), 
'■' rii n 

so the first requirement is met. To control the edge lengths, first remark that 

1{B) = min(6o - ao, 1{B')) L{B) = max(6o - oq, L{B')) 

l{Bij) = min(gi - qi^i,l{B[,j)) L{Bij) = max{qi - qi^i, L{Blj)). 

By the induction hypothesis, the bounds on rij and the fact that L(B') < L{B) we get: 

C{d-l,6) 
^y^ij) ^ i/id-i) ^^^ ) 

On the other hand, reasoning as in the proof of the base case, 

1 77 ■ 

{q^-q^^l)\B'\<-^^l{B) 
n 

KB)<^{bo-ao)\B'\ 



so 



n- 2r)~2 

{q^-q^-l)<{bo-ao)5~^-<L{B) 



Therefore L{Bij) < C {d, 5)n^^' '^ L{B) . The proof of the lower bound on l{Bij) follows the 
same lines and is omitted. This concludes the induction step, and the lemma is proved. D 

Corollary 2.7. Under the assumptions (H1)-(H2)-(H3), for all A C Mi{W^), it holds 
\mi\ni —-K log Z^Qlsii A) > — inf{/(r/); 77 G mt( A), rj <^ Lebesgue}. 

Proof. Let A C A^i(M^) be a Borel set and let rj G int(y4) be absolutely continuous with 
respect to Lebesgue with density h and such that I(rj) < +cx). For all n G N*, define, for 
some sequence e„ converging to 0, 

??„ = (1 - en)l^n + EnAn, 

where dunix) = ^ inm{h{x); n)lr_„.„id(x) dx and d\n{x) = ,^\j lr_„.„id(x) dx, where the 
normalizing constant C„ — > 1, when n — > +00. It holds 

lijln) = (1 - enfK^n) + 2e„(l - £„) / / ^{x, y) dun{x)dXn{y) + ell{Xn). 

Choose En converging to sufficiently fast so that the last two terms above converge to 
when n — )■ 00. Using the monotone convergence theorem it is not difficult to check that 
lit^n) — ^ liv) when n — t- 00 and that r]n converges to rj for the weak topology. 
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Let r > be such that B(ri,2r) C A ; for ah n large enough, B{r]n,r) C B{ri,2r) C A. 
Since rjn satisfies the assumptions of Proposition 2.5, we conclude that for n large enough 

liminf — ^ log ZnQn{A) > liminf — ^ log ZNQNiB{r]n, r)) > -I{rin). 

Letting n — ;► oo and optimizing over {tj £ A,r] <^ Lebesgue} gives the conclusion. D 

End of the proof of Theorem 1.1. The properties of 1^, and the upper bound in point (2) 
are already known. The lower bound of point (2) is given by Corollary 2.7. 

To prove point (3), let A C Mi(W^) be some Borel set and take /U S int(A). According to 
Assumption (H4), there exists a sequence of absolutely continuous probability measures Un 
converging weakly to /x and such that /(z/„) -^ I{^J'), when n — >■ oo. For all n large enough, 
Un ^ A so applying Corollary 2.7, we conclude that liminfAr_!.oo ]^ log.^Ar(5Ar(^) > 
—I{vn). Letting n — )• cxd and then optimizing over /x S int(^) we arrive at 

liminf — ^ log ZnQn{A) > - inf{/(/x); /i G int(^)}. 

According to point (2) of Proposition 2.1, we conclude that Qm obeys the full LDP. D 

2.5. Proof of the almost-sure convergence. Let us establish the last part of Theo- 
rem 1.1. First note that since /^ has compact sublevel sets and is bounded from below, 
/j, attains is infimum, so /min is not empty. For an arbitrary fixed real e > 0, consider the 
complement of the e-neighborhood of /mm for the Fortet-Mourier distance: 

Ae := {IminTe ■= {f^ (^ Ml : d{fl, /„,;„) > e}. 

Since / is lower semi-continuous, Cg := inf^gyi^ /(/x) > 0, thus P(/iAr G A^) < exp(— A^^Cg), 
by the upper bound of the full large deviation principle. By the first Borel-Cantelli lemma, 
it follows that almost surely, limAr_j.oo d{^j\f , Imin) = 0. 

2.6. Sufficient conditions for (H4). The following proposition gives several sufficient 
conditions under which assumption (H4) holds true. Even if some of these conditions are 
quite general. It is an open problem to find an even more general and natural condition. 
One may possibly find some inspiration in [I'lol I]. 

Proposition 2.8 (Sufficient conditions for (H4)). Let V : W^ ^ R and W :R'^ xW^ ^ 
(— cxD,-|-oo] be symmetric, finite on W^ x M.'^ \ {{x,x);x G M.'^} and such that (H2) and 
(H3) hold true. Assumption (H4) holds in each of the following cases: 

(1) W is finite and continuous on W^ x W^. 

(2) For all x G W^, the function y i— )■ W{x, y) is super harmonic, i.e. W satisfies 

W{x,y)> / / W{x,z)dz, Vr>0, 

\B[y,r)\ JB(y,r) 

where \B(y,r)\ denotes the Lebesgue measure of the ball of center y and radius r. 

(3) The function W is such that W{x + a,y + a) = W{x, y) for all x,y,a G M'^ and 
the function J defined by 



J(/i) = J J W{x, y) dfi{x)dfi{y) (2.4) 

is convex on the set of compactly supported probability measures. 

Proof. Let fj, G 7Wi(M'^) be such that /(/u) < oo. First note that, under the assumptions 
(H2) - (H3), the condition /(/i) < oo easily implies that / |F|d/i < -|-oo and // \W\ dfi"^ < 
-hoo. Moreover, it follows follows from (H3) that W{x, y) > (eq - l/2){V{x) + V{y)) + c, 
for some finite constant c. Therefore W is bounded from below on every compact, and so 
the definition (2.4) of J(;u) makes sense if fj, is compactly supported. 
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For all R > 0, let us define fiji as the normalized restriction of fi to [—R; R]"^. Using the 
dominated convergence theorem, it is not hard to see that /i/j converges weakly to fj, and 
that I[^r) — >■ /(/i) when R — )■ +00. To regularize /i/j, we consider /x/j^e = Law(X/j + ef/), 
e < 1, where Xr is distributed according to ^r and U is uniformly distributed on the 
Euclidean unit ball Bi of W^. It is clear that hr^s has a density with respect to Lebesgue. 
Moreover, ///j^^ — )■ /x/j, when e — )■ 0. Indeed, if / : M*^ — >■ M is continuous, it is bounded on 
[-R] R]'^ + Bi, and it follows that j f d^R^e = W{Xr + eU)\ -^ E[/(Xh)], when e ^ 0. 
This applies in particular to f = V . Now let us show in each cases that J[^R^e) converges 
to J{fj,R), when e goes to 0. Let us write 

J(^R,,) = E[W{Xr + eU, Yr + eV)], 

where Xr, Yr, U, V are independent and such that Yr = Xr and V = U. 

(1) If VF is finite and continuous on M x M , then using the boundedness of W on 
{[-R; Rf + Bi) X {[-R; Rf + Bi), it follows that J(/Xi?,£) ^ J{hr) when e ^ 0. 

(2) If W is super harmonic, then We{x, y) := Ej/y [W(a; + eU, y + eF)] < W{x, y) for 
all j;, y. Moreover, it follows from the continuity of W outside the diagonal that, 
for all X / y, We{x, y) — )■ W[x, y) when e — >■ 0. Since /(/x) < +00, ^ does not have 
atoms and so the diagonal is of measure for ^^. It follows from the dominated 
convergence theorem that J[fiR^^) — )■ J[^r) as e — )■ 0. 

(3) Denoting by /ij^ the law of Xr + x, we see that jjlr^s = IE(/[/i/j^]. Therefore, the 
convexity of J yields to 



J(/Ufl,.)<Ec/[J(^f)]=E, 



jj W{x + eU,y + eU) dfiR{x)d^iR{y) 



JifJ-R) 



where the last equality comes from the property W{x + o, y + a) = W{x,y). On 

the other hand, Fatou's lemma implies that liminfe_j.o J^f^R^e) > J{^r). Therefore 

J{f^R,e) -^ JifJ-n), when £ goes to 0. 

We conclude from the above discussion that for any 6 > 0, it is possible to choose R 

sufficiently large and then e sufficiently small so that (i(///j^£,/x) < 5 and [/(//R^g)— /(/x)| < 6. 

This completes the proof. D 

3. Tools from Potential Theory 

In this section we recall results from Potential Theory, that will prove useful when we 
discuss the proof of Theorem 1.2 and Example 1.3. There are many textbooks on Potential 
Theory, with different point of views; our main source is [Lan72] , where the Riesz case is 
well-developed. 

In this section, and unless otherwise stated, we set ka := k/^^ and we take T^(a;, y) := 
ka{x — y), < a < d, d> 1. We denote respectively by M.i C Aioo C A^+ C M.± the sets 
of probability measures, of positive measures integrating ka{-)l\.\yi, of positive measures, 
and of signed measures on M . 

3.1. Potentials and interaction energy. We benefit from the constant sign of the Riesz 
kernel: k^ > 0, contrary to the Coulomb kernel in dimension d = 2 and its logarithm. 
Following [Lan72, p. 58], the potential of /i G A^+ is the function C/^ : W^ — ^ [0, 00] defined 
for every x G M*^ by 

UHix) := fw{x, y) dfi{y) = f k^{x - y) dfi{y). (3.1) 



Note that Ul^{x) = 00 if ;U has a Dirac mass at point x. By using the Fubini theorem, 
for every fi £ A4+, we have C/^ < 00 Lebesgue almost everywhere iff /x G Moo- This 
explains actually the condition < a < d taken in the Riesz potential, which is related 
to polar coordinates {dx = r'^~^drdad). In fact if /x G A^oo then C/^ is a locally Lebesgue 
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integrable function. Moreover, as Schwartz distributions, we have U^ = k^* ^ and, with 
the notations of (1.6), 

-cAall^ = {-cAaka) * /X = /X. 

The interaction energy is the quadratic functional Jq : Ai+ M> [0, cxd] defined by 

Ja(/i) := JJw{x, y) dfi{x)dfi{y) = ju^ d^i. 

Note that Ja{^i) = oo if /i has a Dirac mass, and in particular JaifJ-N) = oo. 

Lemma 3.1 (Positivity and convexity on 7W+). 

• For every /x G M.+ we have JaifJ') > with equality ifffi = 0; 

• Ja '■ -M.^ I-)- [0, oo] is strictly convex: for every fj.,^ £ £+ with fJ- ^ v, we have 

Vt G (0, 1), Jaitfi + (1 - t)u) < tJaifi) + (1 - t)Ja{u); 

• £a,+ '■= {/U S Ai+ : JaifJ-) < oo} is a convex cone. 

We recall that in classical Harmonic Analysis, a function i^ : M x M — )• M is called a 
positive definite kernel when ^^^^XiK(xi,Xj)xj > for every n > 1 and every x G C". 
If this holds only when xi + • • • + x^ = 0, the kernel is said to be weakly positive definite. 
The famous Bochner theorem states that a kernel is positive definite if and only if it is 
the Fourier transform of a finite Borel measure. The famous Schoenberg theorem states 
for every / : M+ — >■ M+, the kernel (x, y) i— )■ /(|x — y\ ) is positive definite on M for every 
d > 1 if and only if / is the Laplace transform of a finite Borel measure on M-|- . The famous 
Bernstein theorem states that if / : M — )• M is continuous and C°°((0,oo)) then / is the 
Laplace transform of a finite Borel measure on M+ if and only if / is completely monotone: 
(— l)"'/^") > for every n > 0. For all these notions, we refer to [BCR84, Kol05]. 

The proof of Lemma 3.1 is short and self-contained. It relies on the fact that the 
convexity of the functional is equivalent to the fact that T^ is a weakly positive definite 
kernel, which is typically the case when VF is a mixture of shifted Gaussian kernels, which 
are the most useful weakly positive definite kernels. For example this works if for some 
arbitrary measurable a, /3 : M — t- M and Borel measure rj, and every x,y G M"', 

POO / 9 2 \ 

Wix, y) = w{\x -y\)= (e-° ^l^"^! + /3(t)) drjit). 

The shift /3 can be < 0, which allows non positive definite kernels such as the logarith- 
mic kernel (note that the Riesz kernel is positive definite). The method is used for the 
logarithmic kernel in [BAG97, Proof of Property 2.1 (4)] with the following mixture: 

1 /■°° 1 / k-y|^ _i\ , 

log = — e 2t — e 2t ] dt. 

This kernel has a sign change and a double singularity near zero and infinity, which can be 
circumvented by using a cutoff. Alternatively, one may proceed by regularization and use 
the Bernstein theorem with the completely monotone function f{t) = (e -|- t)~^, /3,e > 0, 
and then the Schoenberg theorem, see e.g. [Mat97]. For instance, for the logarithmic 
kernel, the following representation is used in [HPOO, Chapter 5]: 

log ,' = f f ^ , - J-) ,t. 

e + \x-y\ Jo \e + l + \x-y\ 1 + tJ 

Finally, let us mention that for the Riesz kernel, yet another short proof of Lemma 3.1, 
based on the formula ka = cka/2 * ^q/2) can be found in [Lau72, Theorem 1.15 p. 79]. 

Proof of Lemma 3.1. Set (5 := d — a. We start from the identity, 

/•oo 

r(l + a)=c^+"/ f'e-'^dt, oO, a>-l. 
Jo 
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Taking c =\x — y\ and 1 + q = 13/2, we get, for every x, y G M"^, 

K{x -y) = J^ fit) e-'\^-y\ dt, where /(t) := ^;^^. 
Now for every n G A^_|_ such that Ja(^) < oo, 

Expressing the Gaussian kernel as the Fourier transform of a Gaussian kernel, we get, by 
writing e*^^"^'*"^ = e*^^'"'^e~*^^''"^ and using the Fubini theorem. 



e -^l^-^l dfi{x)dn{y) = {ATTty^/^ 



A^-y^eiiH' dw] dfi{x)dfi{y) 



^i{x,w) 



dfj,{x] 



e 4t 



rM^ dw. 



Kw{n) 



Now Kw is clearly convex since for every /i, z^ G 7Wi(M ) and every t G (0, 1), 



tKUlJ-) + (1 - t)i^«,(z^) - K^itfi + (1 - t)zy) 



K^(^ - 1^) 



i^i^'^) d{fi - u){x) 



> 0. 



t(l-i) 

It follows then that J^ is also convex as a conic combination of convex function. Let 
us establish now the strict convexity of Jq. Let us suppose that /x, z^ G A4i(M'^) with 
Jailj) < oo and Ja{i^) < oo and tJaifJ-) + (1 — t)Ja{i^) = Ja{tfi + (1 — t)iy) for some 
iG (0,1). Then 

tJo,{fl) + (1 - t) Ja(t/) - Jgjtn + (1 - t)t^) 
t(l - t) 

Arguing as before, we find 



0. 



{Anty^^ [ Ky,{f,-u)e-i->\'dw 



JaifJ- - v) 



{) = J^{^i-u)= I f{t) {ATTty^'' I Ky,{ix-v)e-ttM^ dw dt. 
Jo V Jwi 

Hence, for every i > (a single t > suffices in what follows) 

/ K^(^-i/)e-sl'"l'(iu' = 0. 

Thus, the Fourier transform of ^ — z^ vanishes almost everywhere, and therefore li, = v. 
Finally, £a,+ is clearly a cone, and its convexity comes from the convexity of Jq. D 

Following [LMn72, p. 62], for any // = /i+ — /i_ G M.± such that ^± G A^oo we have 
Ua^ < oo Lebesgue almost everywhere, and we may define for Lebesgue almost every x 

U^ix) := U^+{x) - Uy{x) G (-oo,+oo). 

Following [Lan72, p. 77], for every // = //+— ^_ G A^± such that /i± G M.00 and 

/ U^+ d/i- < 00 and / U^ dfi+ < 00, 

we may define Jaifj) £ (—00, +00] as (thanks to the Fubini theorem) 

J„(^) := jui^ dfi = Ju^+ df,+ + Ju^- dfi. - Ju^+ dfi^ - Ju^- dfi+. 

More generally, for every iJ,i,fi2 G -A^± such that fJ.i±,fi2± G A^oo and 

Ju^'^dfl2:f < 00 



FIRST ORDER GLOBAL ASYMPTOTICS FOR CALOGERO-SUTHERLAND GASES 19 

we may define JaifJ-i, fJ'2) £ (— C)0,+oo] by 

Ja(/il, M2) := Ju^i' dfl2 = JU!^'+ d/X2+ + JU!^'- dfl2- " J U^i'+ d/i2- " J K'- dfl2+. 

Following [Lan72, p. 78], since ka is symmetric, then the reciprocity law holds: 

Ja{fJ-l,fJ'2) = Ja{lJ'2,fJ-l) 1-6. U^^ dH2 = U^^ dfii. 

Let £a be the set of elements of M.± for which Ja makes sense and is finite. As pointed out 
by N.S. Landkof [Lan72] in his preface, a very nice idea going back to H. Cartan consists 
in seeing Ja as a Hilbert structure on £„■ This idea is simply captured by the following 
lemma, which is the analogue of Lemma 3.1 for signed measures of finite energy. 

Lemma 3.2 (Properties of {£a-,Ja))- 

• Ja is lower semi-continuous on Ea for the vague topology (i.e. the weak topology 
with respect to continuous functions vanishing at infinity); 

• £a is a vector space and (^i,//2) ^^ JaifJ-i, 1^-2) defines a scalar product on Sa- 
in particular for every fj, G £a, we have JaifJ-) = JailJ-, m) > with equality iff fi = 0; and 
moreover, Ja : £a ^ (— C)0, 00) is strictly convex: for every 1^,1/ £ £a with ^i i^ y, 

Vy.^/nn ^^M + (1 - t)Ja{y) - Jgjtfl + (1 - t)u) 

VtG(0,l), t(l-t) = Ja(/X - I/) > 0. 

Proof. The lower semi-continuity for the vague convergence follows from the fact that 
ka > 0, see [Lan72, p. 78]. The vector space nature ol £a is immediate from its definition. 
The bilinearity of (/xi,//2) 1— )• J«(/^i,/U2) is immediate. By reasoning as in the proof of 
Lemma 3.1, we get Jq(^, ;^) > for every fj, S £a, with equality iff /i = 0. D 

Following [Lan72, Theorem 1.18 and 1.19 p. 90], for this pre-Hilbertian topology, it can 
be shown that £a,+ is complete while £a is not complete if a > 1, and that Ja is not 
continuous for the vague topology. 

3.2. Capacity and "approximately/quasi every^vhere". The notion of capacity is 
central in Potential Theory. We just need basic facts on zero-capacity sets. Once more we 
follow the presentation of Landkof ([Lan72, Chapter ILl]), to which we refer for additional 
details, references and proofs. 

For any compact set K consider the minimization problem 

WaiK) =inf{J„(i^);z^ G Mi n £:Q,,supp(z^) C K}. 

The boundedness of K implies that Wa{K) £ (0,cxd]. Its inverse Ca{K) is called the 
capacity of the compact set K. The capacity of K is zero if and only if there is no measure 
of finite energy supported in K. 

On general sets on can define an "inner capacity" and an "outer capacity" by 

Ca{A) = sup{Ca{K),K C A,K compact}, 

Ca{A) = inf{C„(0),^ C 0,0 open}. 

It can be shown (see [Lan72, Theorem 2.8]) that if A is a Borel set, these two quantities 
coincide — A is said to be "capacitable" and the common value is called the capacity of 
A. 

A property P{x) is said to hold "approximately everywhere" if the set A = {x, P{x) is false} 
has zero inner capacity, and "quasi-everywhere" if it has zero outer capacity. For many 
"reasonable" P{x), the set A is Borel and the two notions coincide. The following re- 
sult ([Lan72, Theorem 2.1 & 2.2]) shows that, for such "reasonable" properties, "quasi- 
everywhere" means "z/-almost surely, for all measure i' of finite energy". 
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Theorem 3.3. A Borel set A has zero capacity if and only if, for any measure v of finite 
energy, i^(j4) = 0. 

In 'particular, if Ca{A) > 0, A has a positive inner capacity, and there exists a compact 
K C A and a probability measure v of finite energy such that supp(z^) C K. 

3.3. The Gauss averaging principle. In the classical Coulombian case (a = 2), we will 
need the following result, known as Gauss' averaging principle. In W^, for all r > 0, let 
ar be the surface measure on the sphere dB{0, r); its total mass is a^r'^'^ where a^ is the 
surface of the unit sphere. 

Theorem 3.4 (Gauss' Averaging Principle). In M.'^, 

1 /■ 1 J ^:a^ "i-f \x\ < r 

^^^^^, JoBiO,r) ^r^^^'(^) = \^ if\x\> r. 

This result can be found in [Hcl09], Lemma 1.6.1 p. 21. 

4. Proof of the properties of the minimizing measure 

The proof of Theorem 1.2 is decomposed in two steps. We begin by proving the exis- 
tence, uniqueness and the support properties of /i* in Section 4.1. The characterization of 
Hi, is proved in Section 4.2 

Recall that, for a probability measure /U, we have defined 

Hl-i) = -Jail^) + J Vdfl. 

In this section we consider the following minimization problem: 

V:mi{I{fi),fieMi}. (4.1) 

4.1. Existence, uniqueness and compactness of the support. The existence of a 
minimizer for V is clear since we have already seen that / has compact level sets. 
Since I{fi) < oo implies that fj, £ £a and / Vdfj, < oo, the problem V is equivalent to 

Va ■■ inf |/(;u),/i G Ml nSa such that V G -L^(/i)| (4.2) 

in that they have the same values and the same minimizers. Let us call p the common 
value. 

Suppose fj, and u are two measures in Mi H £a such that V G L"'^(/^) H L^{u). Let 
V' : t G [0, 1] h^ [0, cxd) by 

^{t):=I{{l-t)fi + tu) (4.3) 

= i^Jaii'^ - *)/" + t^) + (1 - *) / ^dfi + t / Vdu. 

By Lemma 3.2, ijj is strictly convex ii h y^ v. If ^u and z^ minimize /, then they are in 
A^i n £a, so ■0 is well-defined, and since ■0(0) = /(/i) = /(z^) = V'(l)) A* must be equal to 
1^. Therefore the minimizer //j, is unique. 

Let us now prove that //* has compact support. This result also holds in dimension 2 
with the logarithmic potential, see [ST97, Theorem 1.3, p. 27]. To this end, let us define, 
for any compact K, a new minimization problem: 

Vk :inf {/(^),/i G Mi n <Sa,supp(/i) C K}, 

and let px be the value of Vk- 

Lemma 4.1. Let K be a compact set, and suppose that V{x) > 2p-|-3 when x ^ K, where 
p is the common value of V , Va (defined by (4.1) and (4.2)J. Then the problems V and 
Vk fire equivalent: their values p and px o,re equal, the minimizer exists and is the same. 
In particular, the minimizer /x^ of the original problem V satisfies supp(/i*) C K. 
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Proof. Suppose // is such that I(^) < p+ 1. We will prove that, if n{K) < 1, we can find 
a UK, supported in K such that I[^ik) < I{lj)- This clearly implies that the two values 
Pk and p coincide. Since we know that the minimizer //j, of the original problem exists, 
this also proves that it must be supported in K. 

Let us now construct ^k as the renormalized restriction of /i to K. First, remark that 
n{K) cannot be zero, since 

p+l>I{lJi)>{l-^l{K)){2p + 2,). 

Therefore we can define 

Since by assumption /u(i^) < 1, we may similarly define ^ik"- The measure ji is the convex 
combination 

// = fJ.{K)iJ,K + (1 - ijl{K))iik<'. 
The positivity of V , W, and the choice of K imply that 

I{fl) = haifl) + KK) J VdflK + (1 - KK)) J VdflK^ 

> ^fi{KfJa{fiK) + KKf J VdfiK + (1 - KK)){2p + 3), 
since JaifJ^K") and the interaction energy Jail^Kil^K'^) are both non negative. Therefore 

I{fi) > KKflif^K) + (1 - KK))i2p + 3). 
Assume that I{hk) ^ -^(/")- Then 

I(^)(l-Mi^)2)>(l-^(i^))(2p + 3). 
Using the fact that I{fj,) < p + 1, and dividing by 1 — ij,{K), we get 

2{p + 1) > (p + 1)(1 + fi{K)) > 2p + 3, 
a contradiction. Therefore I{f^K) < I{lAi ^^"^ the proof is complete. D 

4.2. A criterion of optimality. In this section we prove the items (5) and (6) of The- 
orem 1.2. The corresponding result in dimension 2 for the logarithmic potential can be 
found in [ , Theorem 3.3, p. 44], and its proof carries over quite smoothly. 

Proof of item (5) of Theorem 1.2. We already know that /i^ has compact support. The 
first step is to show that ;U^ satisfies (1.7) and (1.8). Let ^ = /i^, and let i/hehw ^ Mir\£a 
such that V G L^(z^). Recall the function ^ from (4.3): 

Since /x^. minimizes /, V'(0"'') must be non-negative. Since Ja is quadratic we get 

= V dii^ + t V d{u - fi^) + -(^J„(/x*) +t^Ja{v - jii,) + 2tJa{lJ.*,i'- li* 
Therefore 

V''(0+) = V d{u - /i*) + JaifJ'*, V - /^*) 

= j{V + U^*)dv-C^ 

= J{v + c/^* - a) du. 
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This must be non negative, so 

Vi/ G TWi n £a, f{V + U^* - C^)du > 0. (4.4) 

Since this holds for ah v, V + C/^* is greater than C^ quasi everywhere. Indeed, let 
A = {x, V{x) + U^*[x) < C^}. Since V + C/^* is measurable this is a Borel set. Suppose 
by contradiction that its capacity is strictly positive. By Proposition 3.3, there exist a 
compact set K <Z A and a measure v with finite energy supported in K. For this measure 
/y + Ul^*dv < C^, which contradicts (4.4). This proves (f.7). 

Let us prove (1.8). Suppose ^(x) + C/^* (x) > C* for some a; G supp(/ij,). Sincey + ?7^* 
is lower semi-continuous, we can find a neighborhood U of x, and an r] > such that 

VxGZ^, V{x) + U^*{x)>C^ + 7]. 

Therefore 

fv + U^i^dfi^ > (a + r/)/i.(Z^) + / {V + U^*)di^^. 



Since V + U^* > C^ quasi everywhere, and /i^ has finite energy, this holds //^ almost 
surely, so 

c\ = Jv + u^*dfi^ > a + ??/u^(z^). 

This is impossible since ii^{U) > 0, by definition of the support. Therefore (1.8) holds. D 

Proof of item (6) of Theorem 1.2. Let v G iSQ,nA^i(M ) such that V G L^(z/) and consider 
again the function '(/;(t) = /((I— t)/u-|-tz/), t G [0, 1]. According to Lemma 3.1, this function 
is convex. Therefore, 

I{u) = ^{l) > V'(O) + ^'(0+) = /(^) + v'(o+). 

To prove that fj, minimizes la, it is thus enough to prove that '0'(O'^) > 0, for all choice of 
v. According to a calculation made above. 



^'(0+) = jv + u^du- lv + u^ d^l 



According to (1.10), V + C/^ > C quasi everywhere. The set where V + U^ < C is Borel, 
so according to Theorem 3.3, it has i' measure 0. Therefore, 



V + U^diy>C. 
On the other hand, (1.9) implies that 

fv + U^dfi= f V + U^dfi = C. 

J Jsupp/^ 

This proves that ^''(0"'') > and completes the proof. D 

Let us sketch an alternative argument for item (6) of Theorem 1.2, based on the so 
called Principle of Domination. This result holds only for a < 2. 

Theorem 4.2 (Principle of domination). Suppose a < 2. Let fi and v he two probability 
measures in 8a, o.nd c a real constant. If the inequality 

C/^(x) < U^ix) + c 

holds fi-almost surely, then it holds for all x G M'^. 

Proof. In the Coulomb case a = 2, [Lan72, Theorem 1.27, p. 110] applies, since U^ is 
super-harmonic. If a < 2, the potential C/" is a-superharmonic, so we can apply [L;\ii72, 
Theorem 1.29, p. 115] and get the result. D 
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Alternative proof of item (6) of Theorem 1.2 in the case a <2. Suppose that // satisfies 
the hypotheses (1.9), (1-10) of the theorem. Let us prove first that C = C*. By (1.9), 
U^ + V = C /x-almost surely, and by (1.7), If^* + V > C^, /x- almost surely. Therefore 
U^ < U^* + C — C^,, /x-almost surely. According to Theorem 4.2, this inequality holds 
everywhere: 

Vx G R'^, C/^(x) < U^* {x) + C-C^. (4.5) 

In particular it holds at infinity where the potentials go to zero (this follows from the 
compactness of the support of /x^, and /x). Therefore C > C^,. Reversing the roles of fx and 
/i^, we find that C must be equal to the "modified Robin constant" C^. 

Coming back to (4.5), we see that U^ < U^* everywhere, and similarly U^* < U^. 
Therefore J{iJ,,fj,) = J{fi,fj,-),) = J(/i*), so 

J(/x - /x*, // - /x*) = J(^, n) - J{iJ,, /x*) + J(/x*, /x*) - J(/x, /x^) 

= 0, 

and n = n-^. D 

4.3. Radial external fields in the Coulomb case — Example 1.3. For the sake of 
completeness, let us finally give a proof of the result mentioned in Example 1.3. Recall 
that V is supposed to be radially symmetric: there exists v : M+ — )■ M such that 

V{x) =v{\x\). 

In the following we suppose that v is C^. Standard computations show that 

AV{x) = ^v'{\x\)+v"{\x\) = -^w'ir), 

\x\ j.a L 

where w{r) = r v'{r). Suppose additionally that AV > 0, i.e. that w is increasing. 

In this case it is natural to look for a radially symmetric equilibrium measure. Guided 
by the results of [' ' ] , we also suppose that /x has a density with respect to Lebesgue 
measure. Let M : ]R_|_ — )■ M+ be compactly supported and let dfj, = M{r)dardr. Suppose 
that the total mass 

roo 

Ud / M{r)r'^~^dr 
Jo 

equals 1. Let us compute the potential of ^. Using the Gauss' averaging principle (Theo- 
rem 3.4), it holds 

U''{x)= ff M{r)W{x,y)dar{y)dr 

M{r) f ^^dar{y)dr 

JdB{0,r) \x - y\ 

= — -T^ / MMr'^-^dr + f3ad / M(r)rdr. (4.6) 

\xf ^ h J\x\ 

To continue, we make the additional assumption that /x satisfies AF = /x on supp /x. In 
other words, we fix a radius i?o and define 

M{r) = C{Ro) (^t^'(r) + v"{r)^ 1||.|<^„} 
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The constant C{Rq) must be chosen so that /i is a probabihty measure, so 



(4.7) 



Since M{r) = for r > i?Oi we have two cases. If |x| > i?Oi the second integral in (4.6) 
disappears, and 



[/^(x) = /3adC(i?o) 



1 



X\ -10 



Ro 



d-l. 



M(r)r'^~'^dr 



I :d-2 "^(Ro) 



If \x\ < Rq, remarking that the derivative of rv'{r) + {d — 2)v{r) is rM{r) (for r < Rq), 
we get 



C/^(x) 



P^Td /-I"! 



\X\" ~ JO 

C{Ro)(3ad 

I \d-2 



M{r)r'^~^dr + Pad j M{r)rdr 

\x\ 



Ro 



Ro 



w{\x\) + C{Ro)pa4rv'{r) + {d - 2)«(r)]q 



= C{Ro)/3a4\x\v'{\x\) + Rov'{Ro) + {d - 2)v{Ro) - \x\v'{\x\) - {d - 2)v{\x\)) 
= C{Ro)pad{Rov'{Ro) + {d- 2)v{Ro) -{d- 2)v{\x\)) 

To satisfy the requirement of the characterization theorem, the term that depends on |x| 
in the last equation has to cancel exactly with the exterior potential V . This occurs if 
C{Ro)l3ad{d — 2) = 1, which can be rewritten using Equation (4.7) as: 

w{Ro)=p{d-2). 

If this holds, the total potential U^ + y is given by: 

[dhi^o^'iRo) + id- '^HRo)) if kl < ^0. 

The potential of fi is such that the hypotheses (1.9) and (1-10) of Theorem 1.2 (item 6) 
hold: we have found the minimizing measure. 

If the external field V is quadratic, i.e. if v{r) = r^, w{r) = 2r'^, so tq = 0, M{r) = 
2dC(i?o)l{|x|<i?o}' ^liere Ro satisfies w{Ro) = 2R^ = f3{d -2),so Ro = {(3{d - 2))V<i. in 
other words the equilibrium measure is the uniform measure on the ball centered in of 
radius (/3(d-2))Vrf. 
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